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ABSTRACT 

We investigate the degree of spatial correlation among extended structures in the 
LMC and SMC. To this purpose we work with sub-samples characterised by differ- 
ent properties such as age and size, taken from the updated catalogue of Bica et 
al. or gathered in the present work. The structures are classified as star clusters or 
non-clusters (basically, nebular complexes and their stellar associations). The radius 
distribution functions follow power-laws (dN/dR oc R~ a ) with slopes and maximum 
radius (R m ax) that depend on object class (and age). Non-clusters are characterised 
by a w 1.9 and R ma x % 472 pc, while young clusters (age < lOMyr) have a ~ 3.6 and 
Rmax ;$ 15 pc, and old ones (age > 600 Myr) have a w 2.5 and R ma x % 40 pc. Young 
clusters present a high degree of spatial self-correlation and, especially, correlate with 
star-forming structures, which does not occur with the old ones. This is consistent 
with the old clusters having been heavily mixed up, since their ages correspond to 
several LMC and SMC crossing times. On the other hand, with ages corresponding to 
fractions of the respective crossing times, the young clusters still trace most of their 
birthplace structural pattern. Also, small clusters (R < 10 pc), as well as small non- 
clusters (R < 100 pc), are spatially self-correlated, while their large counterparts of 
both classes are not. The above results are consistent with a hierarchical star- formation 
scenario for the LMC and SMC. 
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1 INTRODUCTION 

Star formation in the Milky Way and other galaxies is de- 
scribed as a (mass and size) scale-free, hierarchical process, 
in which turbulent gas forms large-scale structures with a 
mass distribution following a power-law. In essence, such a 
scale-free process leads to a mass and size fractal distribu- 
tion. As a consequence, young stellar groupings are clustered 
according to hierarchical patterns, with the great star com- 
plexes (associated with the ~ 10 7 Mq superclouds) at the 
largest scales and the OB associations and subgrou ps, small 
loose groups, clusters and cluster subclumps (e.g. lEfremovl 
1 19951 ) at the smallest. 

In several galaxies, the interstellar gas appears to follow 
a fractal structure ranging from the subpc (« the current 
resolution limit) to the kpc scales; if star formation occurs 
preferentially at the dens est regions, stars should form fol- 
lowing such patterns (e.g. lElmegreen fc Elmegreenll200ll and 
references therein). In this context, star clusters, formed at 
the core (i.e. the densest regions) of giant molecular clouds, 
can be taken as the unavoidable s tar formation pro duct in a 
hierarchically structured gas fe.g. lElmegreenll2006l ). A simi- 
lar picture, in which star clusters are present in dense cores, 



emerges from numerical si mulations that follow in time th e 
collapse of gas clouds (e.g. IWalsh. Bourke fe Myers 2006), 
which also occurs when effe cts of radi ative feedback and 
magnetic fields are included <|Bate!l2009h . 

In a hierarchical scenario, the turbulent gas forms large- 
scale structures (clusters and loose groups) with a mass 
distribution following a power-law of negative slope, i.e. 
dN/dM oc M _/3 , with j3 ~ 2, consistent with th e mass dis- 
tribu tion functions measured in several galaxies (Elmcgrccn 
2008(1 . 

Recent studies came up with robust evidence in- 
dicating that star-forming regions are indeed hierarchi- 
cally structured, for insta nce in the nearby spiral ga laxies 
M33 dBastian etaHl2007h. M51 (jBastian et alj|2005bh . and 
NGC628 (|Elmegreen et al.l I2006T). the Local Group dwarf 
irregular galax y NGC 6822 llKarampelas et al. 20091). the 
Gala ctic disk (Ide la Fuente Marcos fc dc la Fuentc M arcos! 
20091) . and the Gould Belt i|Elias, Alfaro fc Cabrera-Cand 
2002). 



Given the relative proximity, the Magellanic Clouds are 
an excellent environment to investi g ate t he above issues. 
For instance, lEfremov fc Elmegreenl l|l998l ) found that the 
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Figure 1. Top: LMC non-cluster (Table [TJ structures are clumpier (left) than the clusters (right). Bottom: When cluster (red circles) 
and non-cluster (black) angular sizes are shown, hierarchical structuring appears to occur in these typical LMC (left) and SMC (right) 
fields. By far, most of the objects with a large angular size in the bottom panels are non-clusters 



average age difference between pairs of LMC clusters in- 
creases as a function of their distance, which implies hierar- 
chical star formation coupled to evolutionary effects. The 
angular correlation of LMC stellar populations for sepa- 
rations between 2' (~ 30 pc) and 40' (~ 550 pc) also im- 
plies large-scale hierarchical st ructure in current star forma- 
tion (|Harris fc Zaritskvl [l999h . The character of the LMC 
HI structure as a function of scale, the filamentary and 
patchy structures of the high- and low-emission regions re- 
spectively, suggest that most of the ISM is fractal, pre- 
sumably the result of pervasive turbulence, self-gravity, 



and self-similar stirri ng |Elmegreen. Kim fc Stavelev-Smithl 
l200ll ). More recently. iBastian et al I (120091 ) found a highly 
substructured and rapidly evolving distribution in the LMC 
stars. They suggest that all of the original structure is erased 
in ~ 175 Myr (approximately the LMC crossing time), with 
small-scale structures mixing first. Similar conclusions apply 
to the SMC, in which stars appear to have formed with a 
high degree of (fractal) sub-structure, possibly imprinted by 
the turbulent nature of the parent gas; these structures are 
subsequently erased by random motions in the galactic po- 
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tential on a time-scale of a crossin g time through the galaxy 
([Gieles, Bastian fc Ercoranollioosl ). 

In this paper we investigate the degree of spatial corre- 
lation among the different kinds of LMC and SM C extended 
struct ures listed in the updated catalogue of iBica et all 
(2008), together with its relation to star formation. We also 
study properties of their size distribution functions. Only 
two wide-apart age ranges are used for spatial correlation 
purposes: (i) very young objects (not older than ~ 20Myr, 
and probably younger than ~ 10 Myr) , which encompass 
clusters related to nebular emission and associations related 
or not to emissi on, as class i fied a nd catalogued from sky 
survey plates by IBica et ail l|2008h and references therein, 
and (ii) old clusters (older than ~ 600 Myr). According to 
our definition, the dynamical age of the very young clus- 
ters is lower (Sect. [SJ) than the crossing time (of the host 
galaxy), while for the old ones it corresponds to several 
crossing times, which is important for interpreting the spa- 
tial correlation in different time periods. Clusters within the 
wide age range m 20 — 600 Myr are not used in the spatial 
correlation analysis (Sect. Q)). 

Only the Magellanic Clouds have so far such a deep, 
homogeneous information on star clusters, associations and 
nebulae. Exceptions are some neighbouring dwarf galaxies 
that have been surveyed and are (i) featureless (Ursa Mi- 
nor), (ii) contain a few globular clusters (Fornax) or, (Hi) 
star-forming events lik e in the Clouds (e.g. NGC 6822 - 
iKarampelas et aill2009l ). 

This paper is organised as follows. In Sect.[2]we briefly 
discuss the updated Magellanic System catalogue. In Sect. [3] 
we describe the selection criteria for star clusters older than 
the Hyades. In Sect. [4] we discuss the size (and mass, for 
star clusters) distribution functions of the different classes 
of objects. In Sect. [5] we examine the spatial correlation of 
the different structures by means of two-point correlation 
functions. Concluding remarks are given in Sect. [6] 



2 THE UPDATED MC CATALOGUE 

Pr operties of t he up dated MC catalogue are fully discussed 
in IBica et aU (|2008l ). We recall here the basic statistical 
properties. Taking the LMC, SMC and the Bridge together, 
the updated catalogue contains, respectively, 3740 classi- 
cal star clusters, 3326 associations, 1445 emission nebulae, 
and 794 HI shells and supersh ells. With t he rec ent addi- 
tions and cross-identifications, IBica et all (|2008 | ) c o ntain s 
abo ut 12% more o bject s than those in IBica et al.l (|l999l ) 
and lBica fc Dutral l|2000l ) together. 

Especially in view of the spatial correlation analysis 
(Sect. in this paper we restrict the object selection to 
the LMC and SMC, not including Bridge or extended Wing 
structures. A census of the LMC and SMC extended struc- 
tures is provided in Table [T] separated according to ob- 
ject class and including the probable age range. We note 
that, based on similarities observed in the size distributions 
(Sect. fi|, in the present paper we include the AC and NC 
classes (relatively young objects) into the cluster classifi- 
cation, thus res ulting in a higher number of such objects 
than quoted in IBica et ahl ((2008). Besides the latter two 
classes, the cluster classification also contains the C, CN 
and CA classes. As non-clusters (structures mostly associ- 



ated with star formation environments) we take the A, AN, 
NA, DAN and DNC classes (see Table □ notes) . The SNR 
and HI shells are not used because they are object classes 
apart and their size distribution functions are significantly 
different from those of the clusters and non-clusters (Fig. [4]). 

Figure Q] (top panels) shows the angular distribution of 
the 4455 LMC clusters (left) and 2587 non-clusters (right). 
Both kinds of s tructures trace w ell-known LMC (and SMC) 
structures (e.g. IBica et ail 120081 and references therein) . It 
is also clear that the non-clusters appear to present a high 
degree of spatial correlation, with most of them tightly 
clumped together. This applies as well to the clusters, but to 
a lesser degree, because young and old clusters present sig- 
nificantly different levels of spatial self-correlation, the latter 
being essentially non-correlated (Sect.[SJ). 

When the angular sizes are considered (bottom pan- 
els), we see that most structures are arranged according to 
complex patterns, with sub-structures located inside larger 
ones. 



3 OLD STAR CLUSTERS 

The identification, characterisation and spatial distribu- 
tion of old star clusters in the Clo uds has been a m a- 
jor c o ncern throu g hout decades (e.g. | Hodgd [l960l; [Hod 
Il982l ; iBriickl Il975l ; Ivan den Bergbl 1 198 ll ; IBica et al l ll99 
By old or red star clusters we mean those older than the 
Hyadefl (« 630 Myr), or intermediate age clusters (IACs) 
up to classical globular cluster ages. We adop t ed th e def- 
inition of o ld star cluster by Ijanes fc Phelpsl (Il994l ) and 
iFriell i|l995l ). Magellanic Cloud clusters about this age ap- 
pear to show dynamical l y evolved surface density profiles 
(e.g. iMackev fc Gilmorel l2003al ; iMackev fc Gilmorel 12004 
ICarvalho et al.l 120081 ). Since the cluster age distribution 
function drops significant ly with age (see, e.g., Fig. 2 of 
Ide Griis fc Goodwinl [20091 for the Magellanic Clouds), the 
presently adopted old cluster definition encompasses a sta- 
tistically more significant sub-sample (Sect. [372|) than what 
would result for, e.g. clusters older than 1 Gyr. 

Clusters are expected to mix up by random motions 
under the galactic potential on a time-scale of a cross- 



ing time that, for the SMC, is of the order of 75 Myr 
|Gieles. Bastian fc Ercolanol [2008), and about twice that 
value for the LMC0. Thus, the old clusters as defined above 



1 On (blue) sky surveys, both a Hyades-age and a much older 
cluster would appear to consist of a considerable number of stars 
of about the same magnitude. Visually, they would look pretty 
much the same. Young clusters, in contrast, are dominated by 
just a few very bright stars, essentially those at the top of the 
main-sequence turnoff. 

2 As a caveat we note that these dynamical timescales are 
essentially based on random cluster orbits. However, there is 
kinematical evidence suggesting that the LMC cluster system 
rotates as a flattened disk, but the disk geometry and sys- 
temic vel ocity appear to be different for young and old clus- 
ters (e.g. Freeman, Illingworth & Ocmler 1983; Schommcr et al. 
Il992l; iGrocholski et al.ll2006h . Indeed, some studies of the inter- 
mediate age and old populations have found that the velocity 
dispersion increases with age (e.g. [ Hughes. Wood fc Reid|[T99ll ; 
ISchommer et al .1 19921: 1 Graff et al.l2000h . In any case, the dynami- 
cal timescales may be longer than those used in the present paper. 
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Table 1. LMC and LMC Extended Object Properties 









LMC 




SMC 




LMCH 


-SMC 




Class 


Acre 


N 


Rmax 


a 


N 


Rmax 


a 


K 


-L< 1 1 ) CI X. 


Oi 


Comments 




fMvr) 




(pc) 






(pc) 












C 


Any 


2268 


38 


3.53 ±0.22 


456 


30 


3.01 ±0.24 


2724 


38 


3.60 ±0.23 


Ordinary cluster 


CN 


< 10 


81 


25 


3.85 ± 0.62 


9 


8 


2.76 ± 0.58 


90 


25 


3.73 ± 0.67 


Cluster in nebula 


CA 


5-20 


738 


21 


3.19 ± 0.59 


110 


11 


2.59 ± 0.89 


848 


21 


3.10 ±0.57 


Cluster similar to assoc. 


AC 


10 - 30 


1185 


32 


4.42 ±0.24 


60 


11 


2.44 ±0.17 


1245 


32 


4.11 ±0.21 


Assoc. similar to cluster 


NC 


<5 


183 


16 


3.52 ±0.30 


72 


8 


4.15 ±0.66 


255 


16 


3.62 ±0.26 


Nebula w/prob. emb. cluster 


Clusters 




4455 


38 


3.29 ±0.22 


707 


30 


3.04 ±0.23 


5162 


38 


3.20 ±0.19 


C±CN+CA+NC+AC 


A 


< 30 


1476 


171 


2.10 ±0.15 


130 


292 


2.23 ±0.23 


1606 


292 


2.21 ±0.14 


Ordinary association 


AN 


<io 


217 


262 


1.80 ±0.11 


39 


62 


2.07 ±0.30 


256 


265 


1.70 ±0.11 


Association w/nebular traces 


NA 


<5 


817 


472 


1.75 ±0.06 


169 


283 


2.08 ±0.17 


986 


472 


1.73 ±0.07 


Nebula w/embedded assoc. 


DANt+DNCt 




77 


400 


1.15 ±0.09 


33 


144 


1.02 ±0.16 


110 


400 


1.05 ±0.09 


Decoupled structures 


Non-clusters 




2587 


472 


1.94 ±0.06 


371 


288 


2.04 ±0.13 


2958 


472 


1.89 ±0.06 


A+AN+NA+DAN+DCN 


SNR 




52 


78 


0.86 ±0.10 


22 


38 


0.50 ±0.57 


74 


78 


0.85 ±0.13 


Supernova remnants 


HI shells 




124 


472 


3.43 ±0.40 


545 


482 


2.88 ±0.20 


794 


477 


2.82 ±0.05 


HI shells and supershells 



Table Notes. N is the number of objects; a is the power-law slope (<j>(R) = dN/dR oc R~ a ) fitted to the large radii range (Sect. |3J; 
Rmax is the maximum radius measured in each class, (f ) - small cluster or association in large nebula. DCN and DAN: the nebular and 
stellar component of the objects can be distinguished on Sky Survey plates. 



have ages that correspond to several crossing times of the 
respective galaxy, and any memory of the clumpy structures 
where they were born should have been erased. In this sense, 
they can be used as probes of the long-term behaviour of the 
cluster spatial correlation (Sect.0. 



3.1 Short cluster history in the Clouds: towards 
taking the census of the total population? 

iKronl (| 19561 ) and iLindsav 1953 discovered luminous and in- 
ter mediate luminosity cluste rs in the SMC usi ng plate mate- 
rial. |HodgeI&^igM dHzl) and lBriickl (| l975l ) discovered in - 
termediate and low luminosity clusters, while iHodgei (|l986l ) 
discovered even f ainter ones by means of 4m telescope plates. 
iBica fc Schmitq(ll995l ) discovered low luminosity clusters on 
sky survey plates, while IPietrzvnski et all (|l998l ) discovered 
low lumino si ty clu sters by means of CCD imaging. 

IHodgei |l96d ) identified 35 luminous old clusters by 
mean s of non-calibrated CMDs , of which 11 were discov- 
eries. Ishaplev fc Lindsay! (|l963l ) and iLvnga, fc Westerlundl 
i|l963t ) discovered most of the luminous and intermediate- 
luminosity clusters in th e LMC, the latter work being dedi- 
cated to the outer parts. Irlodge fe Sextonj |l966j ) disco vered 
intermediate-luminosity clusters, while Irlodgel ll 19881) low- 
lumin osity ones with 4m-telescope plates. Olszewski et al.l 
ll 19881) discover e d low- luminosity clusters in the outer parts. 
iKontizas et~aT1 (|l990T) discovered additional intermediate 
and low- luminosity clusters, while [Bica et al. (1999) discov- 
ered a l arge number of low lumin osity clusters on sky survey 
plates. IPietrzvnski et all (| 19991 ) discovered low-luminosity 
clus ters in the L MC w ith CCD observations. 

IBica et ah! (|2008l ) and references therein have cross- 



identified these catalogues and a number of other studies, 
and is particularly suitable as a starting point for a deeper 
new survey such as the Visible and Infrared Survey Tele- 
scope for Astronomy (VISTA jf]. It also provides a mean to 
properly acknowledge previous discoveries and to unambigu- 
ousl y establish new cluste r findings. 

ISantiago et alj (Il998l ) serendipitously detected two faint 
clusters in an LMC bar field using HST. The clusters have 
masses comparable to those of Galactic open clusters and 
ages in the range 200 — 500 Myr. The clusters are extremely 
faint on DSS and XDSS images, which suggests that the 
Clouds might harbour an important open cluster counter- 
part population. Besides being a powerful tool to explore 
probable red brighter and intermediate-luminosity star clus- 
ters in the Clouds (Table [2}, VISTA will be essential also to 
detect such a possible population of open cluster counter- 
parts, and estimate their age distribution. 

Based on the broad-band UBVR ph otometry of 
iHunter et all |2003l ). Ide Griis fc Anders! l|2006j ) derived ab- 
solute values of age and mass for a sample of LMC star 
clusters to study the cluster formation rate, their charac- 
teristic disruption time-scale and the cluster mass function 
in different mass ranges. The same sample was used for 
further investigation of the cluster formation rate an d the 
disruption time-scale by IParmentier fc de Gr iis (2008). The 
same method was applied to a sample of SMC clusters by 
Ide Griis fc Goodwin! (|2008l ) to study the infant mortality. 
Our approach in this paper differs in several ways, since we 
intend to build statistically significant samples of clusters 

3 http:/ /www. eso.org/gen-fac/pubs/messenger/archive/no.127- 
mar07 / arnaboldi . pdf 
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(as well as associations and emission nebulae) characterised 
by very different age ranges. 



3.2 Construction of the present sample of old 
clusters 

We compiled ages from the literature later than 1988, as 
determined from CMDs. There are 85 and 202 old clusters 
with ages derived from CMDs for the SMC and LMC re- 
spectively, and they are provided in Table Columns 1 
to 8 of this tab le contain the sa me information as the gen- 
eral catalogue (|Bica et al J 120081 ). We now introduced addi- 
tional columns that provide the age determination method 
(Col. 9), log(Age) (Col. 10), and the relevant references for 
the age (Col. 11). Note that several references are compila- 
tions themselves, so more refer ences are therein. 

We employed observed dRafelski fc Zaritskvl 120051 ). 
reddening-corrected l|Hunter et al.ll2003l ) integrated colours, 
and SW B types to identif y old clusters (typically SWB IVB 
or later. iBica et al.lll996T ). for clusters that still lack CMD 
ages. We also included results from integrated spectroscopy 
|Ahumadaet all 120021 ). By inspection of DSS and XDSS 
images we excluded clusters with apparent contamination 
by relatively bright stars, concerning integrated colours and 
spectra. We found 41 and 117 old clusters in the SMC and 
LMC respectively fro m integrated c olour s (Tabl e 2j). 

Finally, following iBruckl ()l975l ) and iBriickl (|1976l ). we 
examined blue and red DSS and XDSS images, and ESO 
film sky survey plates to identify red clusters. It is remark- 
able how the red SMC clusters bv lBriickl l)l97rJ ) — his types 
Tl and T2 — have been confirmed as old clusters by means 
of deep CMDs. Briick disposed of U plates to help the clas- 
sification. We dispose of blue and red plates, where it was 
basically possible to recognise clusters with brighter stars 
red from the RGB or bluer MS stars. Also, blue clusters 
have as rule more irregular angular distributions. Most clus- 
ters that we examined by this simple method are in the 
outer part s of th e LMC. By means of integrated colours, 
IBica et ail (|l996j ) found that the outer LMC disk appears 
to be essentially composed of old clusters. Our goal here is 
to provide a sample of probable red clusters suitable to cor- 
relation function tests, and to isolate that sample for CMD 
studies in view of VISTA and other large telescopes. There 
are 19 and 203 clusters, respectively in the SMC and LMC, 
that are probably old (red) from our plate inspections. 

Table [2] also includes rather populous clusters that have 
CMDs, integrated colours or plate diagnostics pointing to a 
blue-red transition cluster, that occurs around 500 Myr. The 
sudden or perhaps rather smooth integrated colour change 
is expected from the so-called AGB and RGB phase tran- 
sitions (e.g. iMucciarelli et al.ll2006l and references therein). 
The present sample is a new one for such purposes. To min- 
imise ambiguous age determinations, we have not included 
the blue-red transition clusters in our spatial correlation 
study, although they are, in principle, a lso old enou gh for 



dynamical purposes. Also, we note that IBica et all (| 19961 ) 
decontaminated the clusters containing atypical bright stars 

4 Given the large number of star clusters (667), Table [2] is avail- 
able only in electronic format. Here we provide an excerpt, for 
illustrative purposes. 
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Figure 2. Angular distribution of the LMC (left panels) and 
SMC (right) old star clusters with age obtained by means of 
CMDs, integrated colours and plate inspection. 



supe rimposed. We exami ned all clusters showing red colo urs 
from lHunter et al.l (^OOot ) and lRafelski fc Zaritskvl l|2005h on 
sky survey plates, and excluded those that appeared to be 
dominated by one or a few bright stars. We may have ex- 
cluded some faint intrinsic old clusters with one or a couple 
bright AGB stars, but such sta rs are rare even in populo us 
Magellanic Cloud clusters (e.g. lAaronson fc Mouldlll982h . 

In summary, there are 522 star clusters in the LMC 
that can be currently considered as old as, or older than the 
Hyades. The SMC contains 145 such cases. Considering the 
LMC and SMC together, the total sample of old/red clus- 
ters corresponds to a fraction of ~ 13% of the cluster-like 
structures (Table [5J. Details on the angular distribution of 
this sub-sample are shown in Fig. [2] The old CMD sam- 
ple shows a well-defined LMC b ar, while The SM C prob- 
ably shows a thick edge-on disk (iBica et al.ll2008 j and ref- 
erences therein). Red integrated colours complement these 
samples mostly for fainter clusters. The LMC plate sam- 
ple corresponds essentially to the outer disk. We emphasise 
that the present sharp inner border of the old sample is an 
artifact, but not the outer ring structure, as can be seen 
in the most rece nt census of clusters and related objects 
( Bic a et al . 2008). The outer LMC disk ring is a real feature, 
probably produced as a consequence of the last LMC/SMC 
encou nter that took place « 200 Myr ago (|Bekki fc Chibal 
12007 ) . This struct u re is p resent in the u niform plate surve y 
bv lKontizas et all (|l990h and in that bv lBica et ail l|l999h . 
The m agnitude- limited integr ated photometry of LMC clus- 
ters bv lBica fc Schmitg ll 19951 ) also showed this structure for 
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Table 2. Old SMC and LMC clusters inferred from different methods - Excerpt 



Designations 


a[J2000] 


<5[J2000] Class a 


b 


PA 


Classification 


Method 


log(Age) 


Ref 




(Jims) 


r") o 


(') 


(°) 






(yr) 




(1) 


(2) 


(3) (4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


SMC Star Clusters 



AM-3,ES028SC4 


23:48:59 


-72:56:43 


C 


0.90 


0.90 - 


Old IAC 


CMD 


9.74 


R3 


L1,ES028SC8 


0:03:54 


-73:28:19 


C 


4.60 


4.60 - 


Globular Cluster 


CMD 


9.95 


R12 
















CMD 


9.88 


R26 


L2 


0:12:55 


-73:29:15 


C 


1.20 


1.20 - 




PLA 


Red 


R55 


L3,ES028SC13 


0:18:25 


-74:19:07 


C 


1.00 


1.00 - 




PLA 


Red 


R36 


K1,L4,ES028SC15 


0:21:27 


-73:44:55 


C 


2.20 


2.20 - 




CMD 


9.49 


R18 


BOLOGNA A 


0:21:31 


-71:56:07 


C 


0.80 


0.80 - 


sup 47 Tucanae 


CMD 


IAC 


R56 


L5.ES028SC16 


0:22:40 


-75:04:29 


C 


1.10 


1.10 - 


Old IAC 


CMD 


9.61 


R18 



LMC Star Clusters 



NGC1466,SL1,LW1,ES054SC16,KMHK1 


3:44:33 


-71:40:17 


C 


3.50 


3.50 


Globular Cluster 


CMD 


10.17 


R41 


SL2,LW2,KMHK2 


4:24:09 
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1.60 


1.60 




PLA 


Red 


R55 


KMHK3 


4:29:34 
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0.80 
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4:29:36 
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30 
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SL4,LW4,KMHK7 


4:32:38 


-72:20:27 
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CMD 
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R6 


KMHK6 


4:32:48 


-71:27:30 


C 


0.60 


0.55 


80 


PLA 


Red 


R55 



Table Notes. Cols. 5 and 6: semima jor axes a and b: Col. 7: position a ngle; Col. 9 : old age method. C ol. 10: log(Age) when available, or 
age cl ass. R eferences (Col. 11) : R3: lDa Costal jl999); R 5: iBica et al.l fll996h ; R6: iGeisler et al.l | |l997h - relative ag es; R12: ICrowl et al.l 
d200lh ; R18: |Piatti et alj [1200511 ; R26: lGlatt et al.l IMPOST); R36:lBruckl dl975ft, iBruckl jl976T ); R41: |Piatti et al.l <200a) ; R55: present paper 
- red (old) cluster by plate inspection; R56: iBellazzini, Pancino &: Eerrarol 1120051 ) . 



the oldest age group. In the present study, essentially all 
known red clusters in the outer LMC are included in that 
locus. The geometries of the sub-samples were established 
by each survey, but have apparently not affected the cor- 
relation functions, as shown by the tests with different old 
cluster subsamples (Sect.0. 



In thei r surveys, Bica fc Schmittl (|l995h . | Bica et al.l 

l|l999l ). andU ica fc Dutral (|200d ) employed film copies of 
the ESO Schmidt telescope Red Survey and the UK Schmidt 
telescope SERC-J (blue band) survey in Australi43- The red 
plates trace emission nebulae by means of Ha. The limiting 
magnitudes are R=21.5 and Bj=22.5, respectively. Thus the 
detection limit of stars in clusters is very deep, especially in 
J. However, only the new generation of CCD surveys will 
permit to quantify the completeness of those samples com- 
piled or discovered by our group, and certainly to explore 
an as yet undetected population of fainter objects. 



4 SIZE DISTRIBUTION FUNCTIONS: 
STRUCTURAL HIERARCHY 

The spatial distribution of interstellar gas follows a frac- 
tal structure ranging over many scales, from the subpar- 
sec at the smallest to the cluster and star complexes at 
the largest. This suggests that, if stars are formed mostly 



http: / /www.roe.ac.uk/ifa/wfau/ukstu/platclib.html#UKSTmc 



in the denses t regions, they should also f orm in fractal 
patterns (e.g. lElmegreen fc Elmegreenl 1200 ll and references 
therein). Indeed, the power-law nature of the size distribu- 
tion functi on has been observed in Gala ctic giant molecu- 
lar clouds (Elmcgrccn & Falgaronc 1996) and cloud clumps 
i|Williams. Blitz fc Starklll995l ). The end result of this pro- 
cess is that the cluster size distribution should follow a 
declining power-law with size, a behaviour that has been 
observed in several ga l axies (e.g. lE lmegree n fc SalzerJ Il999l ; 
lElmegreen et al.ll200ll ; iBastian et al.ll2005al )~ 

Our first approach in the investigation of the hierarchi- 
cal structures in the Clouds is by means of the size (i.e. ab- 
solute radius, R c \) distribution functions 4>(Rci) ~ dN/dR c \. 
We build (f>(R c i) individually for all classes of objects listed 
in Table[T]based on the apparent major and minor axes given 
in the updated catalogue, together with the Cloud dist ances 
cZlmc ~ 50 kpc and dsivic ~ 60 kpc (e.g. lSchaeferl l2008). Al- 
though the disks of both Clouds are inclined with respect 
to the line of sight, the effect of the distance correction on 
the absolute radius distribution is small, to within the error 
bars (see App. yQ. 

The radius distribution functions for the cluster-like 
structures (Fig. |3J) are characterised by a steep decline 
for R c \ > 4pc, which corresponds to about 0.25'. As we 
discuss in App. [Bj observational incompleteness probably 
accounts for the shape of 0(-R c i) in the small-size range, 
R c i < 4pc. The maximum radius (Rmax) reached by the 
cluster-like structures in the LMC is R ma x ~ 40 pc, and 
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Figure 3. Radius distribution function of the cluster-like struc- 
tures in the LMC (left panels), SMC (middle), and both Clouds 
combined (right). The radial range not affected by incompleteness 
is fitted with the power-law </>(-R) tx R~ a (solid line). 



Rmax ~ 30 pc in the SMC. Power-law fits (<j>(R c i) oc R c i~ a ) 
to the incompleteness-unaffected range are obtained with 
rather steep slopes, a > 3, especially for the (statistically) 
well-defined distributions. The values of R max and a are 
given in Table [T] 

Figure [4] shows (j>(R c i) for the remaining object classes. 
Except for the SNR, the distributions are qualitatively simi- 
lar to those in Fig. [3] with a decline for large radii. However, 
compared to the cluster-like classes, the power-law slopes are 
significantly shallower (a < 2) and the maximum radii are 
~ 10 times as large, reaching R ma x ~ 500 pc in the LMC 
and Rmax ~ 300 pc in the SMC. The incompleteness-related 
turnover for the A, AN, and NA classes occurs at the same 
radii as that in the cluster-like objects. The HI shells, on the 
other hand, have a turnover at R c \ ~ 70 pc (~ 4.4'), which 
might reflect a real effect, not related to completeness. 

Based on similarities of Rmax and a, we define the C, 
CA, CN, NC, and AC classes as cluster-like structures, while 
A, AN, NA, DNC, and DAN as non-clusters. Their compos- 
ite radius distributions are shown in Fig. [5] together with 
the power- law fit. 

As expected from the above discussion, the cluster-like 
slopes for the LMC, SMC, and LMC+SMC distributions 
(q ~ 3) are significantly steeper than the corresponding ones 
derived for the non-clusters (a ~ 2). Figure[5]also shows the 
distributions obtained by adding all the structures, includ- 
ing the SNR and HI shells. While most of the individual 
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Figure 4. Same as Fig.|3]for the non-Cluster structures. Besides 
being structures apart, SNR and HI shells have distributions sig- 
nificantly different from those of the cluster-like and non-clusters. 

features are preserved, the SMC profile, on the other hand, 
now requires two different power-laws to be described. 

The slopes in the radius distribution of the non-clusters 
are consi stent with those measured for H II regions in spiral 
galaxies l|Oev et aljliobi ). 

4.1 Young and old clusters 

We derive the size distribution functions of the young and 
old star cluster population. We take the CN (age < lOMyr) 
and NC (age < 5 Myr) classes (Table [1]) to represent the 
young star clusters. The old (age > 600 Myr) ones were se- 
lected according to the criteria discussed in Sect. [3] 

Significant differences are observed in the size distribu- 
tion functions (Figure especially in the LMC. The distri- 
bution function of the young clusters falls off with radius at 
a steeper rate than the old ones, reaching a maximum size 
(Rmax ~ 15 pc) less than half of that reached by the old ones 
(Rmax ~ 40 pc). In the statistically more significant distri- 
butions of the LMC and SMC combined (right panels), the 
young clusters fall off with the slope a m 3.6, while the 
old ones have a w 2.5. The differences in Rmax and slope 
probably reflect the several 10 8 yr of dynamical evolution 
of the old clusters, a co nsequence of which is an expansion 
of th e outer parts (e.g. iKhalisi. Amaro-Seoane fc Spurzeml 
20071) of the clusters that survive the infant mortality (e.g. 
Goodwin fc Bastian|[2006l ) phasfl 

6 However, recent evidence suggests that, during the infant mor- 
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Figure 5. Same as Fig.|3]for the composite distribution functions 
of the non-clusters (panels a - c) and clusters (d - f). The sub- 
samples of the young and old clusters are in panels (g) - (i). All 
structures together, including SNR and HI shells, are shown in 
the bottom panels. 



The size distribution functions of the young and old 
clusters fall off at a steeper rate (a > 2) than the non- 
clusters (q < 2). 



4.2 A simple mass distribution function 

Hierarchically structured gas is expected to form star clus- 
ters with mass distributed according to a power-law of th e 
form dN/dM cl oc M c r , with » 2 (e.g. lElmegreerJliioih . 

Below we apply a simple method to analytically trans- 
form the cluster radius distribution function into a mass 
distribution. We wish to test if our sample of LMC and 
SMC clusters basically follow the above mass distribu- 
tion. We caution that our approach to the mass distri- 
bution is a simplification, since we do not take into ac- 
count individual mass-to-light (M/L) ratios, which are 
known to vary considerably between young and old clus- 
ters ( e.g. iBica. Arimoto fc Alloinl Il988l : ICharlot fc Bruzuall 
Il99ll ; iLeitherer et alj Il999l).~ However, the presently ex- 
tracted sample l|Bica et all l2()08h is by far (« 87%) dom- 
inated by young clusters and, thus, large variations of the 
M/L ratio are not expected. Besides, instead of computing 
individual masses from integrated luminosity, we use scaling 



relations that apply well to a wide variety (in terms of age, 
mass, and size) of Galactic star clusters to transform one 
kind of distribution function into another. In this sense, we 
expect that the adopted radius to mass transformation is 
representative of the average M/L ratio of the clusters. 

We start by assuming a spherical star cluster with a 
mass radial density profile that can be described by a King- 
like functiorQ (jm{R) = 0"a/o/(1 + (R/R c ) 2 ), where <tmo is 
the surface mass density at the cluster centre and R c is the 
core radius. We also consider that essentially all stars are 
contained within ^ R ^ Rci, where 7? c i is the cluster 
radius. The spatial mass density of such a structure can be 
computed from inversion of AbelPs integral, 



p(R) = -- 
n 



<?M0 

"2~r7 



r d * Vx 2 - R 2 

Thus, the cluster mass can be computed from 
Rci 

Mci » / p(R)Atv R 2 dR = 2-rgmoR 2 x 



1 + (R/Rc 



3/2 



arcsinh(R c \ / R c ) 



1 



1 + (r c /r c1 ) 2 



Galactic star clusters with ages from a few Myr to 
~ lGyr, masses within 50 M < M cl < 7 x 1O 3 M , and 
radii within 2pc < R c i < 20 pc, have the relation be- 
tween R cl and R c well approximated by R c \ ?a 9R C (e.g. 
iBonatto fc Bical l2009ah - Comparable r atios are observed 
in LMC and SMC star clusters (e.g. IMackev fc Gilmorel 
l2003al ; IMackev fc Gilmore!l2003bl ; ICarvalho et al.ll2008h . Un- 
der these assumptions we have M c \ ~ A-komoRc 2 ~ 



densities in the range 3OM pc 



o-jvfo < 6OOAf pc" 



accou nts for the distribution of cluster mass and core ra- 
dius IIBonatto fc Bical l2009bl ). Then, the transformation 
of the radius distribution to mass is given by c/>(M c i) = 
<f)(R c i)/ (0.31a mo Rci), and the average cluster mass-density 
is a declining function of the cluster radius, p(M pc~ 3 ) = 
« ^M^-i. For a given R c 



the radius to 



(4/3)7ri? c ?' 

mass scalings depend only on a mo as M c \ oc a mo 



and 



4>(M C \) oc o M0 that, for different values of cta/o, preserve 
the shape of the mass distribution, only changing the mass 
values. 

We use the above scaling relations to transform the clus- 
ter radius distribution functions (Fig. [5} into mass distribu- 
tions, <f)(M c \)dM c \ = <f)(R c \)dR c \, separately for the LMC 
and SMC, and the combination of both, LMC+SMC. We 
also consider the age ranges < 10 Myr, > 600 Myr, and clus- 
ters of al l ages combined. The a verage value (Galactic star 
clusters - IBonatto fc Bicall2009bl ) of the central mass density, 
a mo ~ 3OOM pc -2 , is used to compute the mass distribu- 
tions (Fig. [6}; fit parameters are given in Table [3] The mass 
distribution of the very young clusters in both Clouds falls 
off at a steeper rate towards large masses than that of the 
old ones, which is c onsistent with a ma ss-dependent disrup- 
tion time-scale (e.g. lLamers et aT]|2005l ). Also, the slopes are 



tality, the star cluster populatio n has been depleted by less than 
» 30 %, both in the the SMC (Ide Griis fc Goodwin! I2OO8I) and 
LMC Jde Griis fc Goodwinll2009fK 



Similar to the function introduced bv lKind lll962t) to describe 
the surface brightness profiles in the central parts of globular 
clusters. 
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Figure 6. Estimated mass distribution functions with a power- 
law (4>(M C \) oc M c i~' 3 ) fitted to the large cluster mass range of the 
LMC (left panels), SMC (middle), and the combined LMC+SMC 
(right) clusters. Age ranges are < lOMyr (top panels), > 600 Myr 
(middle), and all ages (bottom). 



steeper in the LMC than in the SMC. These slopes are con- 
sistent with those of the mass distrib utions of star clusters 
in different galaxies (|Elmegreenll200ct ). 

The mass distributions of the very young clusters are 
characterised by a maximum mass of M max w 1.2xlO 4 M 
(LMC) and M max » 3 x 10 3 M Q (SMC), while for the old 
ones it is M max w (7- 8) x 10 4 M in both Clouds. The lat- 
ter values are considerably higher th an the maximum mas s 
of typical Galactic open clusters (e.g. iPiskunov et al.ll2007l ). 
The decline in the number of clusters with mass below M m i n 
(Table |3j) is probably related to observational incomplete- 
ness in the detection of small clusters (Sect. [4J . For compar- 
ison purposes we also show in Fig. [6] the mass distributions 
for the LMC+SMC clusters, as well as those corresponding 
to clusters of all ages, in which the basic features of the 
individual distributions are preserved. 

Interestingly, the maximum mass of the LM C clusters 
younger than w 30 Myr (Ide Griis fc Andersll2006l). and S MC 
ones younger than « 10 Myr |de Griis fc GoodwirJl2008T ). is 
about 2.4 times higher than that of the very young LMC 
and SMC cluste rs (Table O. For LMC clusters younger 
than « 5.6 Gyr (Ide Griis fc Anders! 120061) and S MC ones 
younger than « lGvr (|de Griis fc Goodwin! 120081 ) . the ra- 
tio increases to « 4. Although the somewhat different 
age ranges, consistency between both sets of M max val- 
ues can be reached with the central mass densities <tmo ~ 
700M Q pc -2 and o mo ~ 12OOM0pc 2 , respectively for the 
very young and old clusters. The somewhat higher values least for a maximum separation £ m ax (which should scale 



of umo in the MCs clusters, with respect to the Galac- 
tic open clusters, is consistent with t he relative cluster 
mass ranges e ncompassed by the LMC (Ide Griis fc Anders! 
l200fj). SMC dde Griis fc Goodwin! liooil ), and Milky Way 
( Piskunov et al. 2007) mass distributions. 

Finally, if we take into account variations of M/L 
with age for the dominan t (in number) young clusters (e.g. 
iBruzual fc Cha riot 2003), the actual mass values for clus- 
ters younger than ps 8 Myr would be > 30% lower than the 
average-(M/L) estimates above, and w 20 — 40% higher for 
those with age within ~ 13 — 30 Myr. G iven that the number 
of clu sters decreases with age (see, e.g. Ide Griis fc Goodwin! 
2009, for the age-distribution of the MC clusters) , our mass 
estimates in each bin of the mass distributions (Fig. [6} may 
be somewhat overestimated. 



5 HIERARCHY ASSOCIATED WITH STAR 
FORMATION 

In a hierarchical scenario, young star clusters are expected 
to preserve some memory of the physical conditions prevail- 
ing in their birthplace. Because of random motions along 
many orbits under the galactic potential, the spatial distri- 
bution of old star clusters, on the other hand, should be 
very little reminiscent of the primordial one. According to 
this scenario, the frequency of young star clusters lying rel- 
atively close to each other - and to star-forming structures 

- should be higher than for the old ones. Based on the 590 
LMC clusters (~ 13% of the pr esent sample size - Table [1} 
catalogued bv lBica et all (Il996t) wit h ages derived from the 
UBV colours bv lGirardi et al.l (jl995l ). lEfremov fc Elmegreenl 
( 1998) found that the average age difference between pairs 
of clusters increases with the separation, which they inter- 
preted as resulting from star formation that is hierarchi- 
cal in space an d time. A similar result - and interpretation 

- was found by Ide la Fuente Marcos fc de la Fuente Marcos! 
(2009) for pairs of open clusters in the Milky Way disk. 

We investigate this point further by means of the degree 
of spatial correlation among groups of objects characterised 
by different age ranges and sizes. We use the young and old 
star clusters defined in Sect. 14.11 

Consider two groups of objects, A and B. For each 
object in A, we compute the angular separation with re- 
spect to all objects in B. After applying the same proce- 
dure to all objects in A, we build the two-point-correlation 
function (2PCF), which measures the fractional number of 
pairs N that lie within a given separation £ and f + d£ , 
2PCF(C> = dN/d^. According to this definition, the 2PCF 
is simply the angular separation distribution function. 

Artificial 2PCFs built with samples of points that emu- 
late both the geometry and object distribution of the LMC 
and SMC are used to check the statistical significance of 
the spatial correlations. In the simulations we randomly se- 
lect the right ascension (a) and declination (S) coordinates 
of a given point within the actual ranges spanned by each 
Cloud (Fig. [TJ and with the same number-frequency as the 
observed ones. 

Irrespective of the adopted geometry, a random distri- 
bution of objects would produce a number of neighbours 



within a given separation £ that increases as N(£) oc £ , at 
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Table 3. Mass distribution properties 



LMC 



SMC 



LMC+SMC 



Age range 
(Myr) 



N M„ 



M„ 



N M„ 



Mmax 



N M min 



Mm, ax 



(M ) (M ) 



(M ) (M ) 



(M ) (M ) 



< 10 201 5.0 x 10 2 3.2 x 10 4 2.24 ±0.14 
>600 435 9.0 x 10 2 7.0 x 10 4 1.74 ±0.11 
All ages 3700 5.0 x 10 2 7.0 x 10 4 1.92 ±0.12 



47 5.0 x 10 2 3.0 x 10 3 1.83 ±0.16 
135 5.0 x 10 2 8.0 x 10 4 1.40 ±0.08 
629 5.0 x 10 2 8.0 x 10 4 1.81 ±0.13 



198 5.0 x 10 2 3.2 x 10 4 2.32 ±0.15 
552 9.0 x 10 2 8.0 x 10 4 1.69 ±0.10 
4271 5.0 x 10 2 8.0 x 10 4 1.87 ±0.10 



Table Notes. N is the number of clusters used to fit the mass range M m i n < M c \ < M„ 
and M ma x computed for trjfo = 3OOM pc~ 2 ; they scale linearly with &mo- 



with the power-law 4>(M C \) oc M a \ P. M m i n 
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Figure 7. Two-point correlation functions for the LMC (left pan- 
els) and SMC (right) extended objects. The simulated 2PCF (dot- 
dashed line) increases linearly with the separation §■ Panels (a) - 
(d): spatial self-correlation for the non-clusters and clusters. Pan- 
els (e) - (f): degree of spatial correlation of the young and old 
(older than the Hyades) clusters with the non-clusters. Panels (g) 
- (h): same as above for the clusters older than 1 Gyr, with the 
age determined from CMDs. 



with the angular size of the simulated field). Thus, the 
2PCFs should increase with £ as dN/d£ oc £. Indeed, the 
2PCFs derived with the simulations (Fig. [7)l present the ex- 
pected dependence with separation for £max < 80' ft! 1200 
(LMC) and t\ max < 40' w 700 pc (SMC). Beyond these val- 
ues both the measured and simulated 2PCFS consistently 
drop, as a consequence of the limited size of the Clouds. 

As a first step we compute the spatial self-correlation 
functions, in which A = B, for the LMC and SMC (Fig. 0. 
Compared to the simulated 2PCFs, the non-clusters (top 



Figure 8. Same as Fig. [7] for the composite distribution func- 
tions of the non-clusters (top panels) and clusters (bottom), but 
differentiating for object size (R). 



panels) present a relatively high degree of spatial self- 
correlation for separations smaller than f < 15' as 220 pc 
(LMC) and £ < 25' w 440 pc (SMC). Young (age < 10 Myr) 
clusters present a high degree of spatial self-correlation, from 
small to large scales (middle panels), £ < 35' ~ 500 pc 
(LMC) and £ < 25' » 440 pc (SMC). Old (age > 600 Myr) 
clusters, on the other hand, have a very low degree of spa- 
tial self-correlation, restricted to separations £ < 0.6' w 9pc 
(LMC) and £ < 1.7' « 30 pc (SMC). Interestingly, the same 
pattern is obtained with the 2PCFs computed for the clus- 
ters older than 1 Gyr, with the age determined from CMDs. 
Some degree of spatial correlation at small separations 
among old clusters is expected, since the Clouds contain 
binary and/or merger star clusters preferentially of compa- 
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rable ages (e.g. Bica et al, 1 19991 ; iDieball. Miiller fc Grebell 
|2002| ; ICarvalho et all l2008h . In summary, young clusters 
have a probability of being clustered together significantly 
higher than old ones, both in the LMC and SMC. The dy- 
namical age of clusters older than 600 Myr corresponds to 
> 4 crossing times in the LMC, and > 8 in the SMC, while 
for the young ones (age < 10 Myr) it is < 0.07 and < 0.13, 
respectively for the LMC and SMC. Given that a single 
crossing time is necessar y to smear out most of the primor - 
dial structural pattern (iGieles. Bastian fc Ercolanol [2008). 
the above conclusion is consistent with the old clusters hav- 
ing been mixed up by the random motions under the galactic 
potential along several 10 8 years, while the young ones still 
trace most of the birthplace pattern. 

Now we test the degree of spatial correlation of the 
young and old clusters with the non-clusters (star-formation 
environments). As expected from the self-correlation analy- 
sis, the young clusters are highly correlated with the non- 
clusters (bottom panels). The old clusters, on the other 
hand, appear to have some spatial correlation with the non- 
clusters only at the very-small scales, £ < 0.2' « 3 pc (LMC) 
and £ < 0.4' » 7pc (SMC). Part of this correlation may be 
due to projection effects on the bar. For larger separations 
the 2PCFs can be accounted for by the random distribution 
of objects. 



5.1 Self-correlation and object size 

Now we examine the spatial self-correlation among clusters 
and non-clusters of different radius ranges. Based on the 
respective radius distribution functions (Fig. 0), we take 
R = 10 pc as the boundary between small and large clusters; 
for the non-clusters we take the boundary at R = 100 pc. 
The derived correlation functions (Fig. [8} indicate that the 
small clusters are more spatially correlated than their large 
counterparts. The same applies to the non-clusters. 

Again, this picture is what should be expected from a 
hierarchical structure. 



5.2 Effective separation 

Finally, we investigate the effective separation of the young 
and old clusters with respect to the non-clusters. We first 
compute the effective separation (£ e //) between a cluster 
and a non-cluster, which we define as the ratio of the angular 
distance (£, converted to the absolute scale) to the radius of 
the non-cluster (Rnc), £eff = £/ 'Rnc- In this way, a cluster 
that is located inside a non-cluster has £ e // < 1. 

Two-point correlation functions built with all the £ e // 
between samples A and B thus provide a measure of the 
clustering among objects in both samples. For comparison 
we use simulated 2PCFs built for object samples (young and 
old clusters and non-clusters) with coordinates selected as 
described in Sect. [5] We now also include absolute radii sep- 
arately for each sample, randomly taken from the observed 
distributions (Fig. [5} ■ 

The resulting 2PCFs are shown in Fig. [9] Consistently 
with the analyses of the previous sections, young clusters in 
both Clouds present a high degree of clustering with the non- 
clusters, especially for small effective separations (£ e // < 
1), but reaching as well high values of £ e //- In all scales, 
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Figure 9. Two-point correlation functions of the effective sepa- 
rations (§ e // = f / Rnc ) between the young (top panels) and old 
(bottom) clusters with respect to the non-clusters, for the LMC 
(left panels) and SMC (right). Simulated 2PCFs are also shown 
(heavy-solid line). £ e // = 1 is indicated by the dashed line. 



the clustering degree of the old clusters with respect to the 
non-clusters, on the other hand, can be accounted for by a 
random distribution of old clusters. 

The above results are consistent with a strong hierar- 
chical structuring of the young star clusters in both Clouds, 
including their time evolution effects. 



6 SUMMARY AND CONCLUSIONS 

In broad lines, when star formation occurs in turbulent gas, 
large-scale structures are expected to be produce d following 
a pow er-law mass distribution (dN/dM tx M~ 2 - lElmegreenl 
120081 ) , an d with hierarchically clustered y oung stellar group- 
ings (e.g. lEfremovlll995l ; lElmegreenlbOOrj ). 

In the present paper we address the above issue by in- 
vestigating the degree of spatial correlation among sets of 
LMC and SMC extended structures, characterised by differ- 
ent properties, and its relation t o star formation. Based on 
the catalogue o f B ica et all (|2008T l , we built sub-samples that 
basically contain star clusters (young and old) and nebular 
complexes (and their stellar associations). The latter struc- 
tures are related to star-forming regions; for simplicity, we 
refer to them as non-clusters. 

In all cases (Figs. 131 oil . the radius distribution func- 
tions follow a power-law (dN/dR oc R~ a ) decline for large 
radii with slopes that depend on object class (and age). Tak- 
ing both Clouds combined, the non-clusters fall-off with a 
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slope a ~ 1.9 and reach sizes of R m ax < 472 pc. Old (age 
> 600 Myr) clusters present the somewhat steeper slope 
a ~ 2.5, while the young (age < 10 Myr) ones have the 
steepest slope a ~ 3.6. The maximum size reached by clus- 
ters is less than w 10% of the non-clusters, with the old 
ones reaching a size « 3x bigger than the young ones. The 
differences in slope and maximum size between the young 
and old clusters can be accounted for by long-term dynam- 
ical effects acting on the clusters. By means of a radius to 
mass scaling (Sect. l4~2"1l . we show that the mass distribution 
of the LMC and SMC clusters follows dN/dM cl cx M A ~ P 
(Fig. 0, with (3 « 2. Within the uncertainties (Table EH, 
this value agrees with the slope expected in a hierarchical 
scenario (Elmegreen 2008). Also, the mass distribution for 
clusters younger than ~ 10 Myr falls off towards large masses 
faster than the clusters older than w 600 Myr. 

According to the two-point correlation functions 
(Sect. the LMC and SMC star clusters younger than 
~ 10 Myr present a very high degree of spatial correla- 
tion among themselves and, especially, with the non-clusters 
(Fig.0. Clusters older than the Hyades (> 600 Myr) on the 
other hand, appear to have been heavily mixed up, proba- 
bly because their ages correspond to several galactic crossing 
times and the strong pertu rbations associated wi th the LMC 
and SMC encounters (e.g. lBekki fe Chibal [20071 ). When the 
analysis is restricted to clusters older than 1 Gyr - with the 
age determined from CMDs, the same conclusions are ob- 
tained. 

Considering two different radius ranges, we show that 
small clusters (R < 10 pc) and non-clusters (R < 100 pc) are 
spatially self-correlated, while the large ones are not (Fig.|SJ). 
Also, young clusters in both Clouds present a very high de- 
gree of spatial clustering with the non-clusters, which does 
not occur with the old ones (Fig. |5J). 

The above results, expressed in terms of the spatial and 
size distribution of extended structures in the LMC and 
SMC, are fully consistent with a hierarchical star-formation 
scenario, in which star complexes are part of a continu- 
ous star-formation hierarchy that follows the gas distribu- 
tion. Similar conclusions drawn from different methods and 
samples of objects have bee n obtained for the LMC (e.g . 
Elmegreen fe Efremovl Il996l; Efremov fe Elmegreen] 1 19981: 
Harris fc Zaritskvlll999l; iLivanou et alj|2006h and the SMC 



fe.g. lLivanou et al.| [2007; Gic les. Bastian fc Ercoland i2008) 



VISTA and other large telescopes will certainly uncover 
a large number of faint clusters in the Magellanic Clouds, 
with masses comparable to the Galactic open clusters. The 
same for embedded clusters. CMDs will provide accurate 
ages for them, as well as for many luminous, intermediate- 
luminosity, and low-luminosity cluste rs already catalog ued. 
Also in this context, the catalogue bv lBica et al.l (|2008T ) will 
be an essential tool, for having gathered and cross-identified 
the small and large structures in the Magellanic Clouds. 
The catalogue is also useful to help establish discoveries. 
The present work has provided as well a sub-catalogue of 
old and probable-old clusters, which can be useful also for 
VISTA studies. 
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APPENDIX A: DISK INCLINATION 

The LMC and SMC disks are inclined with respect to the 
line of sight, and this should introduce some variations in 
the absolute sizes across the disks, as compared to the no- 
inclination approach adopted in Sect. [4] 

To investigate the inclination effe ct on the radius distri - 

^ ]l99Ch 



icr 



bution funct ions we assume 42° (e.g. Kontizas et alj 



and 40° (e.g. lStanimirovic. Stavelev-Smith fc Jonesll2004l ') 
the inclination of the LMC and SMC disks with respect 
to the line of sight. Then, the absolute size of each cluster 
was recomputed for its corrected distance. The inclination- 
corrected radius distribution function (for the combined 
LMC+SMC clusters) is shown in Fig. IA1I in which the un- 
corrected distribution (Fig. [S] panel f ) is also shown for com- 
parison purposes. We conclude that differences are small, 
essentially within the error bars. This can be accounted for 
by the relatively small distance corrections (with respect to 
the adopted Cloud distances), and that corrections affect 
objects both in the near and far sides in the opposite sense. 
On average, when a large number of objects is considered, 
the near and far-side corrections tend to self-compensate. 

Another effect that might introduce variations on abso- 
lute cluster size is the triaxial nature of the Clouds. The 
SMC, for ins tance, may have a line-of-sight depth of 6 
- 12kpc (e.g. ICrowl et al1l200ll ). Thus, depth corrections 
would be of the same order as those related to inclination. 



APPENDIX B: SURFACE BRIGHTNESS 
INCOMPLETENESS EFFECTS 

Since extended structures are the focus of this work, the sur- 
face brightness (SB) incompleteness - which is expected to 
affect the radius distribution functions - should be taken into 
account. Basically, for a given luminosity, a more extended 
object will have on average a lower SB, and thus may not 
be detected by depth-limited surveys. 

We examine this effect by means of a sample of 
10 r artificial star clusters whose luminosity and ra- 
dius distributions are described by 4>(L) dL oc L~ 2 dL 
an d 4>(R) dR oc R~ 3 3 dR, res pectively. As discussed 
in lElmegreen fc Elmegreenl (|200lf ), these analytical func- 
tions describe star clusters and Hll regions. To re- 
produce the input radius and luminosity distributions, 
the simulated radius and luminosity are computed from 
R = Rm/ [l + ni {{R m /R M ) 2/i - I)] 172 ' 3 and L = 
L m / [1 + n%(L m jLu - 1)], where R m , R M , L m , and L M 
are the minimum and maximum radii and luminosities, 
and ni and ri2 are random numbers in the range [0.0, 1.0]. 
The radius and luminosity of a given cluster are indepen- 
dently assigned. This process allows that clusters of the 
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Figure Al. The inclination-corrected radius distribution func- 
tion (filled symbols) is very similar to the uncorrected one (empty 
symbols). 



same size (and mass) - but different ages - may have dif- 
ferent luminosities, as expected from the fading lines as- 
sociated with the stellar evolution. Also, clusters with any 
radius within (R m ,RM) are allowed to have any luminosity 
within (L m , Lm). Then, the SB is computed in the usual 
way, fi — — 2.5 log(L/7r R 2 ) + cnt, in arbitrary units. 

The results are summarised in Fig. IB1I In general, the 
SB distribution among the simulated clusters agrees with 
the expected relation of decreasing SB with cluster radius 
(panel a). The two power-laws that describe the radius and 
luminosity distributions are reflected on the shape of the SB 
distribution (panel b), which first (beginning at the smallest 
and most luminous clusters) increases exponentially towards 
lower SBs, reaches a maximum, and falls off exponentially 
towards the largest and less luminous clusters. Based on 
this distribution we arbitrarily apply cuts for clusters with 
H < 17, 12, and 10 and compute the corresponding radius 
distributions (panel c). Clearly, the SB cuts preserve the 
power-law character of the radius distribution. The major 
effect is a steepening of the slope. lndeed,while the complete 
radius distribution is a power-law of slope a = —3.3, the SB- 
restricted distributions have a = —3.6 ± 0.1, — 5.0 ± 0.1, 
and a = —5.2 ± 0.1, respectively for [i < 17, 12, and 10. 

In summary, SB-related incompleteness affects the ra- 
dius distributions preferentially at the large-clusters tail, 
having little effect on the small clusters. Also, it preserves 
the power-law character of the radius distribution and, due 
to the preferential effect on large clusters, it produces a 
steepening of the distributions. Thus, the decrease in the 
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Figure Bl. Panel (a): model SB distribution with respect to clus- 
ter radius; the density of points (number of simulated clusters) 
roughly increases towards heavier shades of gray, (b): SB distri- 
bution function (fitted with exponentials - solid line) showing the 
arbitrary thresholds (dashed lines), (c): radius distributions cor- 
responding to the SB cuts in (b). 



observed radius distributions towards small clusters (Figs. [3] 
-[SJ cannot be accounted for by SB incompleteness, and ap- 
pears to be linked to an observational effect. In this sense, 
VISTA will be important also to explore the structure of 
small clusters in the Clouds, and to investigate the shape of 
the radius distribution at the small scales. 



